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What entanglement is present in naturally occurring physical systems at thermal equilibrium? 
Most such systems are intractable and it is desirable to study simple but realistic systems which can 
be solved. An example of such a system is the ID infinite-lattice anisotropic XY model. This model 
is exactly solvable using the Jordan- Wigner transform, and it is possible to calculate the two-site 
reduced density matrix for all pairs of sites. Using the two-site density matrix, the entanglement of 
formation between any two sites is calculated for all parameter values and temperatures. We also 
study the entanglement in the transverse Ising model, a special case of the XY model, which exhibits 
a quantum phase transition. It is found that the next-nearest neighbour entanglement (though not 
the nearest-neighbour entanglement) is a maximum at the critical point. Furthermore, we show 
that the critical point in the transverse Ising model corresponds to a transition in the behaviour of 
the entanglement between a single site and the remainder of the lattice. 

PACS numbers: 03.65.Ud, 73.43.Nq, 05.50.+q 



I. INTRODUCTION 

It seems to be a truism in quantum physics that 
strongly entangled systems exhibit complicated be- 
haviour which is difficult to quantify. Two practical 
examples of this 'principle' are the conventional super- 
conductor |^] and the fractional quantum Hall effect 
(FQHE) H . In both cases, for certain parameter regimes, 
the system enters a very interesting entangled state (the 
BCS ground state for the superconductor ||, |5|, and the 
Laughlin ground state for the FQHE |^]). For many years 
these systems resisted attempts to understand them us- 
ing reasoning based on classical methods [[30| . It required 
a major breakthrough, the construction of an insightful 
ground state ansatz, to elucidate the physics of both the 
FQHE and the superconductor. The key feature of both 
systems, which makes it hard to explain them classically, 
appears to be that their ground states arc strongly en- 
tangled. 

Entanglement is a uniquely quantum property of any 
nonlocal superposition-state of two or more quantum sys- 
tems jj], @, H . Such states are typified by the Bell state 
l^ - ) = (| 01) — 1 10)). The many curious features of 
entangled states have motivated considerable research. 
A remarkable consequence of this work is the emerging 
understanding of entanglement as a resource [l^] , like 
energy, which can be used to accomplish interesting phys- 
ical tasks. 

The similarities between entanglement and energy ap- 
pear to be more than just superficial. It turns out to be 
possible to quantify the entanglement present in a given 
quantum state. This allows the development of quan- 
titative high-level principles governing the behaviour of 
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entangled states, independent of their particular physical 
representation. These principles can be seen as analogous 
to the laws of thermodynamics governing the behaviour 
of energy, independent of the specific form in which it 
is given to us. We hope that the quantitative theory 
of entanglement may provide a powerful unifying frame- 
work for the understanding of complex quantum systems. 
This is because, when viewed in terms of their entan- 
glement content, a large number of apparently different 
states turn out to be equivalent. 

This paper is one step in testing the hypothesis 
Jll[ [T^ , fl3[ |T4| ] that the study of complex quantum sys- 
tems may be simplified by first analysing the static and 
dynamic entanglement present in those systems. We will 
attempt to perform such an analysis in a representative 
system chosen from condensed matter physics, specifi- 
cally, the XY model fTq ]. The signature of complexity in 
this system is the occurrence of a quantum phase transi- 
tion. 

Quantum phase transitions (QPT) are a qualitative 
change in the ground state of a quantum many-body 
system as some parameter is varied |l6| , Unlike or- 
dinary phase transitions, which occur at a nonzero tem- 
perature, the fluctuations in a QPT are fully quantum. 
Typically, at the critical point in parameter space where 
a QPT takes place, long-range correlations in the ground 
state also develop. The existence of a QPT in a quan- 
tum many-body system strongly influences the behaviour 
of the system near the critical point, with the develop- 
ment of long-range correlations and a nonzero expecta- 
tion value for an order parameter Jl6| . 

In Jllll it was argued that QPTs are genuinely quantum 
mechanical in the sense that the property responsible for 
the long-range correlations is entanglement. It was also 
argued that the system state is strongly entangled at the 
critical point. It would be desirable, to begin with, to 
show that systems near quantum critical points can be 
simply characterised in terms of their entanglement con- 
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tent. Unfortunately, such a proof seems very difficult. 
We need first to understand the entanglement in such 
systems before proposing a classification scheme based on 
entanglement content. At the moment the most promis- 
ing technique to study entanglement in critical quantum 
systems appears to be the renormalization group, which 
is the standard way to obtain information about systems 
at and near criticality. 

The renormalization group (RG) is based on the no- 
tion that physics at small length scales (and hence higher 
energy scales) should not affect physics at much larger 
length scales. The RG is, in fact, a family of methods 
which can be applied to learn non-perturbative infor- 
mation about strongly interacting systems. The devel- 
opment of the renormalization group (see, for example, 
p8| , p"9| ] for a review) has shown that phase transitions are 
universal in the sense that many properties of the sys- 
tem do not depend on the detailed dynamics of the sys- 
tem under consideration. Instead, using RG techniques, 
it has been shown that phase transitions depend only on 
certain global properties, such as symmetry and dimen- 
sion. We would like to apply the ideas of the RG to 
calculate entanglement quantities in systems exhibiting 
a quantum phase transition. To sec if this is possible, 
it is desirable to first carry out exact calculations in or- 
der to determine if similar universality properties govern 
the entanglement present in such systems. The purpose 
of this paper is therefore to do such calculations for the 
XY model. 

Unfortunately the modern theory of entanglement (see, 
for example, the review issue of Quantum Information 
and Computation |2^, |^, ^3|) is only partially devel- 
oped, and at the present time can only be applied in a 
limited number of scenarios. In these limited scenarios 
well-developed analytic tools exist to quantify the struc- 
ture of entanglement present in a system. Two important 
scenarios are (a) the case of a pure state of a bipartite 
system, that is, a system consisting of only two compo- 
nents; and (b) a mixed state of two spin-i particles. 

For this reason, we focus our investigation on two types 
of calculation for the XY model. The first calculation is 
of the entanglement between a single site in the lattice 
and the rest of the system, for the ground state of the 
model. The second calculation is of the entanglement be- 
tween two sites of the lattice at arbitrary temperatures 
and separations, allowing us to determine whether there 
are truly quantum features present in the two-body cor- 
relations in the system. Thus, although we do not obtain 
an understanding of the three- and more-party entangle- 
ment present in the system, we do calculate significant 
partial information characterising the entanglement. 

The entanglement present in condensed-matter sys- 
tems has been investigated previously by a number of 
authors @, g g ||, pi |7[ |§ [|| g |l[ |§ . It was 
considered by Nielsen |12| who studied the Heisenberg 
model on two sites analytically. An expression for the 
ground state entanglement in the infinite ID Heisenberg 
chain was obtained soon after by Wootters [p5| . Numeri- 



cal calculations of entanglement in the Heisenberg model 
on a small number of sites were carried out by Arnesen 
et al. ]2l| . Arnesen et al. identified parameter regions 
where there is appreciable thermal entanglement, which 
is entanglement present at nonzero temperatures. Re- 
cent studies include the numerical calculation of entan- 
glement in the transverse Ising model on small numbers 
of sites [^8| , and analytic computations of entanglement 
in the XY model on 2 sites Q and 3 sites @. Addi- 
tional studies have been carried out on itinerant fermion 
systems |T^] and other small condensed matter systems 
related to the XY model 



The structure of this paper is as follows. In Sec. II the 
exact solution and calculation of the correlation functions 
for the XY model is outlined using the Jordan- Wigncr 
transform. The thermal ground state properties of this 
system are considered in Sec. Ill, focusing on the special 
case of the transverse Ising model, and the role entangle- 
ment plays in the quantum phase transition in this model. 
Thermal entanglement in the transverse Ising model is 

We conclude in Sec. ^, 



then calculated in Sec. IV 



and 



sketch some possible future research directions. 



II. EXACT SOLUTION OF THE XY MODEL 

In this section we consider the exact solution of the 
XY model on TV sites, which is facilitated by use of the 
Jordan- Wigner transform |3j| . The observables that are 
important for the calculation of the entanglement are 
evaluated in the large-N or thermodynamic limit. The 
two fundamental objects constructed in this study are 
the one- and two-site density matrices. From knowl- 
edge of these matrices it is possible to calculate the one- 
and two-party entanglement occurring in the XY model. 
The solution of the XY model is well known, and the 
procedure outlined in this section to solve it follows the 
standard method |H| |ll| 34, 



The main result in 
this section is the explicit construction of the one- and 
two-party density matrices for the XY model at thermal 
equilibrium. 

The Hamiltonian for the anisotropic XY model on a 
ID lattice with N sites in a transverse field is given by 

il 



H 



N-l 

£ 

3=0 



[(1+7)«J«J +1 + (1-7W + J + 



(1) 



where <r" is the ath Pauli matrix (a = x, y or z) at site j, 7 
is the degree of anisotropy, and A is the inverse strength of 
the external field. We assume cyclic boundary conditions, 
so that the TVth site is identified with the Oth site. The 
standard procedure used to solve Eq. (Q) is to transform 
the spin operators cr° into fermionic operators via the 
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Jordan- Wigncr transform 



where 



d = Y[[-°f]<Ti , 

3=0 
3=0 



(2) 
(3) 



where 



4 = \^ + ^ Vi=\{<r*-i°l). (4) 

It is easy to verify that Cj satisfy the fermionic anticom- 
mutation relations 



{c i ,c]} = S ij , {ci,Cj} = 0. 



(5) 



In terms of the fermionic operators, Eqs. (g)-(^), the 
Hamiltonian Eq. (hi) assumes the quadratic form 



(N-l N-l \ 

Y, c tAi, j c j + -J2 (ctBycj + h.c.) 
i,j=0 i,j=Q J 



N, 
(6) 



■1, A 



57A — ^4_i+i,i , -Bj,i+i 



where A^j 

— Bi + i t i = and all the other -Aj^ and Bi :J 

are zero. The quadratic Hamiltonian Eq. (g) may be 
diagonalised by making a linear transformation of the 
fermionic operators 



N-l 

1 



h q ic\ 1 , 



i=0 
N-l 



4 = E 



(7) 

(8) 



w 9 = v /( 7 Asin^) 2 + (l + Acos^) 2 , (14) 

and g = 2irq/N. 

Now that the Af Y Hamiltonian has been diagonalised 
we can calculate the one- and two-site density matrices. 
Much of the remainder of this paper is concerned with 
the case where the system is at thermal equilibrium at 
temperature T. The density matrix for the XY model 
at thermal equilibrium is given by the canonical ensem- 
ble p = e- fm /Z, where (3 = l/k B T, and Z = tr(e~' 9H ) 
is the partition function. The thermal density matrix is 
diagonal when expressed in terms of the Jordan- Wigncr 
fermionic operators i] q . Our interest lies in calculating 
the quantum correlations present in the system as a func- 
tion of the parameters /?, 7, A. In general this problem 
requires knowledge of all the possible spin correlation 
functions. These correlators are typically very difficult 
to calculate from p as it is diagonal in terms of the r] q s, 
which arc complicated nonlocal functions of the original 
spin operators. Fortunately, the only correlation func- 
tions which we require are the one- and two-point corre- 
lation functions. The evaluation of these functions has 
been carried out previously |?5|, |3?J . 

The one- and two-site density matrices may be con- 
structed from the one- and two-point correlation func- 
tions, using the operator expansion for the density matrix 
of a system of N spin-^ particles in terms of tensor prod- 
ucts of Pauli matrices. For the single-site density matrix 
pi for the first spin — equal, by translational symmetry, 
to the state pi of a single spin at an arbitrary site — the 
operator expansion reads 



pi = tr ;(p) 



(15) 



where q = -N/2, -N/2 + 1, . . . , N/2 - 1 and the g qi 
and h q i can be chosen to be real. By requiring that the 
operators r\ q obey fermionic anticommutation relations, 
and that the Hamiltonian Eq. ([!]) be manifestly diagonal 
when expressed in terms of the fermionic modes r\ q , the 
following two coupled matrix equations must hold 



(A + B)V g = uj q $ q , 



(9) 
(10) 



where the components of the two column vectors $ q and 
ty q are given by 



= 9qi + hqi, 
[®q]i = 9qi - h qi- 



(11) 

(12) 



The quadratic Hamiltonian Eq. (^), when expressed in 
terms of the operators rj q , takes the diagonal form 



(13) 



where trj is the partial trace over all degrees of freedom 
except the spin at site i, o~f are the Pauli matrices acting 
on the site i with the convention a° = Ii, and the coeffi- 
cients q a are real. The coefficients q a are determined by 
the relation 



q a = tr(cr» = (a l a ). 



(16) 



To completely specify the single-site density matrix re- 
quires knowledge of three expectation values (qo = 1 be- 
cause pi must have trace unity). However, because the 
Hamiltonian for the XY model Eq. ([!]) possesses symme- 
tries it is possible to reduce this number to one. First of 
all, the Hamiltonian is real, so that p\ = pi- As the ma- 
trix o~ v is imaginary this means that q2 must be zero. The 
second symmetry that the XY Hamiltonian possesses is 
the global phase flip symmetry 



N-l 



Upf = n o~j. 



(17) 



3=0 
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This symmetry implies that [o z ,pi] = 0, so forcing (73 
to be zero. The single-site density matrix p\ is therefore 
determined solely by q\ . 

For the two-site density matrix, which is the joint state 
of two spins at sites % and j, the operator expansion takes 
the form 



Pij 



E 



3 

a,/3= 



:0 Paf3(7i 



The coefficients are determined by the relation 
tr(ofofpy) = (of of ), 



Pa/3 



(18) 



(19) 



so that if the relevant correlation functions are known it 
is possible to construct the two-site density matrix com- 
pletely. 

The operator expansion Eq. (^) implies that we need 
sixteen correlation functions to construct the two-site 
density matrix. However, as in the case of the single-site 
density matrix, this number can be reduced by appealing 
to the symmetries of the the Hamiltonian. Translational 
invariance of the lattice means that the density matrix 
depends only on the distance r = \j — i\ between the 
spins, that is, p^ = po r . Reflection symmetry about any 
site also means that p^ = pj t . Also, since the Hamilto- 
nian is real, p^ * = pij . Finally, the the global phase flip 
symmetry implies that [o?o|,pjj] = 0. The symmetries 
of the XY model require that the only nonzero coeffi- 
cients in the operator expansion Eq. ( |l8| ) are poo, P03, 
P30, P11, P22, and P33. Furthermore, poo = 1 because the 
density matrix must have trace unity, and P03 = P30- 

In the thermodynamic limit, N — * 00, sums that ap- 
pear in the expectation values are replaced by integrals, 
and the correlation functions for the XY model can be 
reduced to quadratures 15,3^, |37], Bq] . The calculations 
are rather involved, and we merely summarise the results 
here. In thermal equilibrium, for arbitrary 7 and A, the 
transverse magnetisation (o z ) is given by [B5[ 



M = -I 

7T 



dip (1 + A cos ( 



tanh(i/3o;0) 



(20) 



where we abuse notation and write lo^ = LJ q to indi- 
cate the replacement of <p q with the continuous variable 
<j) which results from the thermodynamic limit <j) q — > (f>. 
The two-point correlation functions are given by |37|] 



G—i G-2 


■ • G — r 


Go G_i 


■ ■ G_ r +i 


G r -2 G r _3 


•• G_ x 


Gi Go 


G- r +2 


G2 G\ 


G-r+3 


G r G r —\ • • 


Gi 



(a z a*) = 4(o 2 ) 2 - G r G„ r , 



(21) 

(22) 
(23) 



where 



G r 



d(j) cos(0r)(l + Acos( 

7A r 

d(f> sin(0r) sin(0) 



^tanh(i/3w0) 
UJ4, 

tanh(i/3w0) 

^4, 



(24) 



Summarising, in the thermodynamic limit we may 
write the single-site density matrix pi entirely in terms 
of the transverse magnetisation Eq. (EQ), 



Pi 



(25) 



Similarly, the two-site density matrix po r can be written 
entirely in terms of the correlation functions Eq. (pi]), 
Eq. (p2|), Eq. ( p3| ) and the transverse magnetisation, 

Por = ; • (.^oj 



III. GROUND STATE ENTANGLEMENT FOR 
THE TRANSVERSE ISING AND XY MODELS 

In this section we discuss the quantum correlations oc- 
curring in the ground state of lattice systems undergoing 
a quantum phase transition. We argue that the criti- 
cal point corresponds to the situation where the lattice 
is critically entangled, where, somewhat loosely, we de- 
fine critically entangled to mean that en tangle ment is 

we out- 



present on all length scales. In subsection III A 
line the properties of the ground state of the transverse 
Ising model, which is a sim ple su bclass of the anisotropic 
XY model. In subsection 1IIB the contribution to the 



ground-state correlations from one- and two-party entan- 
glement in the XY model is calculated explicitly in order 
to illustrate the sharp peak in the entan glement at the 
critical point. Finally, in subsection III C we discuss how 
the properties of shared entanglement may be related to 
critical quantum lattice systems. 

In |ll[] it was argued that the physical origin of the 
correlations which occur in systems exhibiting a quantum 
phase transition is quantum entanglement. We reproduce 
the argument of (ll ] here in order that this study be self- 
contained. For concreteness, we restrict our attention to 
a lattice of spin-i particles. 

Suppose the ground state of a quantum lattice system 
were not entangled, that is, it is a product state. Then 
a simple calculation shows that the spin-spin correlation 
function (of of) — (of) (of) is identically zero. Thus, if 
the correlation function is non-zero then the ground state 
must be entangled. Furthermore, we conjecture that 
large values of the correlation function imply a highly 
entangled ground state; it is interesting open problem to 
prove a precise form of this conjecture. 

For general quantum lattice systems the correlation 
function decays exponentially as a function of the sepa- 
ration \i — j\ when the system is far from criticality p6[. 
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When the system is at a critical point, the correlations 
decay only as a polynomial function of the separation. 
At this point a fundamental change in the ground state 
has occurred. 

We believe that when a system approaches a critical 
point the structure of the entanglement in the ground 
state undergoes a transition. Further, we conjecture that 
the nature of this transition is governed by a change in 
the spatial extent of the entanglement. The entangle- 
ment between a single spin and the rest of the lattice 
away from the critical point must be bounded in finite 
regions because the correlations are damped exponen- 
tially. At the critical point correlations develop on all 
length scales, and the physical property responsible for 
these correlations, entanglement, should become present 
at all length scales as well. We believe that a fundamen- 
tal transition in the nature of the entanglement in the 
system occurs at this point; in some sense, at the critical 
point the state is delocalized, compared with the local 
nature of the entanglement away from the critical point. 
If this physical picture is correct, there should be evi- 
dence of entanglement developing on all length scales in 
the one- and two-party entanglement results. 

As described in detail below, the ground state of the 
XY model exhibits the features we have described in 
the previous paragraphs. That is, maximality of the en- 
tanglement at criticality, and evidence that a transition 
in the entanglement structure takes place at the criti- 
cal point. Although much work remains to be done to 
flesh out this physical picture, we believe that further 
research will show that these are generic properties of 
critical quantum systems. 



A. Properties of the transverse Ising model ground 

state 

The ground state of the XY model is very compli- 
cated with many different regimes of behaviour |3^, ^7) . 
For the sake of clarity, we focus most of our discussions 
on the transverse Ising model, which arises as the zero- 
anisotropy limit 7 — * 1 in Eq. (Q). The reason for this 
particular choice is because the transverse Ising model is 
the simplest quantum lattice system to exhibit a quan- 
tum phase transition jlTJ]. The central goal in this sec- 
tion is to illustrate the intimate relationship between the 
entanglement structure of the ground state and the quan- 
tum phase transition. In particular, the calculations for 
the transverse Ising model provide the clearest evidence 
for the conjecture that the critical point corresponds to 
the situation where the lattice is most entangled. 

The Hamiltonian for the transverse Ising model may 
be obtained from the XY model Hamiltonian, Eq. (Q), 
by setting 7 = 1: 

N-l 

a 1>^;<-i •<*;)• (27) 

3=0 



The structure of the transverse Ising model ground state 
changes dramatically as the parameter A is varied. The 
dependence of the ground state on A is quite complicated. 
However, it is possible to investigate the A = and A — > 
00 limits exactly. 

When A approaches zero, the transverse Ising model 
ground state becomes a product of spins pointing in the 
positive z direction, 

|0> A _ tO w.--|T)il T>i+i--" - (28) 

In the A — > 00 limit the ground state again approaches a 
product of spins pointing in the positive x direction, 

|0+> A ^oo « ■ ■ ■ I ->> J -+i - ■ ■ • (29) 

The A — > 00 limit is fundamentally different from the 
A = case because the corresponding ground state is 
doubly degenerate under the global phase flip, Eq. JT?]), 
where 

|(T)a^oo = C/pf|0 + ) a ^oo pa--- 1 Oil Hj+i • • • (30) 

is a second ground state. The A = ground state is 
invariant under the global phase flip. We note that in 
both limits the ground state approaches a product state. 

Using the solutions obtained for the limiting cases of 
A we can qualitatively describe the ground state as A is 
varied. When A is small, the exchange term (TjcrJ +1 may 
be regarded as a perturbation, and perturbation theory 
may be used. In this case the ground state becomes a 
superposition of the unperturbed ground state and low- 
lying excitations in such a way that the small-A ground 
state remains invariant under the global phase flip. 

When A is much greater than one, 1/A is a small pa- 
rameter and perturbation theory may again be used to 
show that the now-degenerate ground states are a super- 
position of the unperturbed ground states |0 + '~) and low- 
lying excitations. The degeneracy of the ground state un- 
der the global phase flip remains for A large. (This degen- 
eracy, along with the invariance of the ground state |0) 
under Upf may be established nonperturbatively [fjTjf .) 

When A = 1 a fundamental transition in the form of 
the ground state occurs. The symmetry under the global 
phase flip breaks at this point and the system develops a 
nonzero magnetisation (o~ x ) 7^ which grows as A is in- 
creased. The magnetisation is the order parameter which 
identifies the existence of a new phase. 

Now that we have outlined the structure of the ground 
state for the transverse Ising model as a function of A we 
have a basic physical picture with which to interpret the 
exact results. 

The calculation of the entanglement between a single 
site and the rest of the lattice requires construction of 
the single-site density matrix for the ground state. While 
the single-site density matrix for the thermal state was 
constructed in Sec. II, there is a distinction between the 
zero-temperature limit of the thermal density matrix and 
the ground state, because of the possible ground-state de- 
generacy. In the following, when referring to the ground 
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state of the system, we suppose the system to be in one 
of the possible degenerate eigenstates |0 + ) or |0~) rather 
than any other linear combination. It does not matter 
which of the two is chosen to be 'the' ground state be- 
cause all the entanglement quantities calculated in this 
paper do not depend on the choice, due to the local sym- 
metry connecting the two states. Therefore, without loss 
of generality, when the system is in the ground state we 
choose the system to be in the eigenstate |0 + ) for A > 1 
and |0) for A < 1. For simplicity we will identify |0 + ) 
with |0) when A is greater than or equal to one. 

The zero temperature state, po, of the XY model may 
be found by taking the limit f3 — * oo of the canonical 
ensemble, 



,-/3H 



Po 



lim ■ 



(31) 



When the ground state is nondegenerate the zero temper- 
ature state is the same as the ground state of the system, 
po = |0)(0|. However, if the ground state is degenerate 
the zero temperature ensemble becomes an equal mixture 
of all the possible ground states. For the transverse Ising 
model the zero temperature state may be written 



Po = ||0+><0+| + ||0-><0-|. 



(32) 



In order to differentiate between the actual ground state 
|0) of the XY model and the zero temperature ensemble 
we refer to po as the thermal ground state. 

In general, the canonical ensemble ^possesses the same 
symmetries as the Hamiltonian Eq. (|lj). This is a simple 
consequence of the identity [U, H] = 0, where U is some 
unitary or antiunitary operator representing the symme- 
try operation. The invariance follows from [U, p] = 0, 
so that U pW = p. In particular, while each individual 
degenerate ground eigenstate may not possess the same 
symmetries as the Hamiltonian, the thermal ground state 
Po has all the same symmetries. 

The quantum phase transition in the transverse Ising 
model separates two different phases, the paramagnetic 
phase where the magnetisation (a x ) is zero, and the 
ferromagnetic phase where the magnetisation becomes 
nonzero. Associated with the development of a nonzero 
value for the order parameter (o~ x ) is the breaking of the 
phase flip symmetry. The symmetry breaking present in 
the ground state |0) is a key feature of the quantum phase 
transition, and is responsible for the development of non- 
zero order parameter (a x ) associated with the ferromag- 
netic phase. (In practice, small external perturbations 
force spontaneous symmetry breaking of the phase flip 
symmetry, and the system will choose one or the other 
ground state, so this order parameter is, in principle, ob- 
servable.) This symmetry breaking cannot occur in the 
thermal ground state. For this reason, we will be most 
interested in properties of |0) rather that pq. For each 
of the degenerate ground eigenstates |0 + ) and |0~) the 
global phase flip symmetry is broken, so the terms which 
were set to zero in the operator expansion Eq. (|l§|), as 



a consequence of the symmetry Eq. (|17|), may become 
nonzero. 

The single-site density matrix p\ for the ground state 
of the Ising model is obtained by taking a partial trace 
over all but one site of |0)(0|. In general, because the 
global phase-flip symmetry may be broken, the operator 
expansion for p\ is only constrained by the reality con- 
dition p\ = p\. Therefore, typically, two parameters are 
required to specify p\ completely, the magnetisation (o~ x ) 
and the transverse magnetisation (cr z ): 



I 



Pi = 



(33) 



It is difficult to calculate the magnetisation (a x ) of the 
ground state explicitly because its expression in terms 
of Jordan- Wigner fermions is nonlocal, but it is possible 
to obtain (a x ) from the large-r limit of the correlation 
yielding 



function {a x a x +r 




A > 1. 



(34) 



The transverse magnetisation (a z ) is given by the integral 
Eq. (^0|) which reduces to an elliptic integral for 7 = 1 
and (3 — » oo, 



1 



do 



1 + A cos < 



v/l + A 2 + 2A< 



(35) 



Armed with knowledge of the appropriate correlation 
functions we can now proceed to the calculation of the 
entanglement in the ground state of the XY and trans- 
verse Ising models. 



B. Ground state entanglement in the transverse 
Ising model 

Given the modern understanding of entanglement as a 
physical resource it makes sense to ask how much entan- 
glement there is in a given multipartite state. In order 
to answer this question the notion of an entanglement 
measure has been developed. A review of work on entan- 
glement measures may be found in [^o[ ^lj ^2], ^3) . 

The study of entanglement measures is far from com- 
pletely developed. There is currently no consensus as to 
the best method to define an entanglement measure for 
all possible multipartite states. There are, however, sit- 
uations where there is an unambiguous way to construct 
suitable measures. It is these situations that we study in 
this paper. 

When a bipartite quantum system AB is in a pure 
state there is an essentially unique measure of the en- 
tanglement between the subsystems A and B given by 
the von Neumann entropy S (7[ |3^, ^o|, El]]. The von 
Neumann entropy is calculated from the reduced density 
matrix pa or ps according to the formula 



S = -tr(p A logpA) = -tr(p B logp B ). 



(36) 
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When either subsystem A or B is a spin-i system, S 
varies from (product state) to S = 1 (maximally en- 
tangled state). For the ground state of the transverse 
Ising model we regard a single site as subsystem A and 
the rest of the lattice as subsystem B. 

When a bipartite system AB is in a mixed state there 
are a number of proposals for measures of the entangle- 
ment in the state, including, the entanglement of forma- 
tion |7], [2^], the distillable entanglement |7], ^2), and the 

. Each of these 
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relative entropy of entanglement 
measures have the property that, for pure states of AB, 
they reduce to the von Neumann entropy. The entangle- 
ment of formation &{A : B) is the best understood of 
the mixed-state entanglement measures. For this reason, 
in this paper, we use the entanglement of formation to 
measure the mixed-state entanglement in the XY model. 

At the current time, there is no simple way to calculate 
the entanglement of formation for mixed states of bipar- 
tite systems AB where the dimension of A or B is three 
and above. However, for the case where both subsystems 
A and B are spin-^ particles there exists a simple formula 
from which the entanglement of formation can be calcu- 
lated . In this case the entanglement of formation is 
given in terms of another entanglement measure, the con- 
currence C |2^, ^6|. The entanglement of formation 
varies monotonically with the concurrence. 

The entanglement S between a single site and the rest 
of the lattice represents the collective contibutions of the 
entanglement between the given site and all other sites 
in the lattice. Unfortunately the single-site entanglement 
does not tell us how the entanglement is shared out. For 
example, 5 = 1 could mean that the site in question is 
maximally entangled with a neighbouring site, or, entan- 
gled with many sites. In the transverse Ising model it 
appears that S is related to the onset of correlations in a 
fairly direct way (see below), and to reflect this we speak 
of S as 'measuring' how entangled the lattice is. 

We should point out that this situation is by no means 
typical. It is quite common for the ground state of 
a condensed matter system to possess strong nearest- 
neighbour entanglement and no long-range correlations 
(see, for example, the models constructed by Affleck, 
Kennedy, Lieb and Tasaki discussed in [[l7| ) . Analysis of 
the entanglement in various AKLT models carried out by 
the authors has shown that, in fact, the single-site entan- 
glement is constant for all parameter values even though 
long-range correlations develop and vanish. The entan- 
glement in these models (and many other condensed mat- 
ter systems) is, in general, not revealed from knowledge 
of the single site density matrix. What is really needed 
- but which has not yet been developed — to study 
these models is an entanglement measure which can take 
account of the way entanglement is shared out. 

At the critical point, A c = 1, of the transverse Ising 
model there is a fundamental transition in the struc- 
ture of the ground state. The correlation function 
(erf)((7^) decays polynomially as a function of 




WTO 



0\ 



separation at this point (the dominant term has exponent 



FIG. 1: Single-site entropy S for the thermal ground state po 
(solid) and the single-site entanglement for the ground state 
|0 + ) (dashed) of the transverse Ising model 



— j) while for all other values of A this decay is exponen- 
tial. Interestingly, one could argue that the correlation 
function itself actually constitutes an entanglement mea- 
sure for pure states as it transforms as a tensor under 
local unitary operations and is zero for product states. 
As argued earlier, the change in the correlation func- 
tion signals a fundamental change in the entanglement 
present in the ground state. This change is reflected in 
the single-site entanglement S for the ground state which 
appears in Fig. [l| The single-site entanglement varies 
from zero at A = 0, where the ground state is a product, 
to a maximum at the critical point A — 1. As the limit 
A — > 00 is approached S also approaches zero because the 
ground state again approaches product form. The single- 
site von Neumann entropy for the thermal ground state 
of the transverse Ising model is also shown in Fig. [l]. Un- 
like the ground state case, the entropy approaches unity 
in the limit A — > 00. This is because the thermal ground 
state approaches an equal mixture of two pure states (the 
cigenstatcs |0 + ) and |0~)) in this limit. The single-site 
entropy is not measuring the entanglement content of the 
thermal ground state in this limit, rather it is measuring 
the degree of mixedness of the thermal ground state. 

It is an intriguing fact that systems with quite dif- 
ferent microscopic dynamics may behave cquivalently at 
criticality. Further, their behaviour depends only on the 
dimension of the system and the symmetry of the or- 
der parameter. The character of this behaviour is cap- 
tured by a small number of universal quantities whose 
behaviour at criticality is completely described in terms 
of a unique single number, a critical exponent. The equiv- 
alence of physically different systems and their simple 
dependence on certain global properties at criticality is 
known as universality. One of the triumphs of twentieth 
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century physics was the development of the renormaliza- 
tion group (RG), which provided an explanation for the 
emergence of universality in critical systems. 

If we are to suppose that S is a universal quantity 
which could be studied via the RG then we should be able 
to find a critical exponent for S. In other words, near the 
critical point we should be able to write something like 



S oc |A — A r 



(37) 



where 7 is the critical exponent for S. Unfortunately, 
this is not possible. As we describe below, the single- 
site entanglement is two-sided, so that two numbers are 
needed to specify S like Eq. (|37|) near the critical point, 
one for each of the two ways of approaching A c = 1. In 
this way we see that the single-site entanglement is not 
a universal quantity. 

The two-sided behaviour of the single-site entangle- 
ment arises because the single-site density matrix de- 
pends on both the magnetisation and the transverse mag- 
netisation. In the region near A < 1 only the transverse 
magnetisation is nonzero and the single-site entropy rises 
linearly. At the critical point the magnetisation becomes 
nonzero and increases as As . This becomes the dominant 
term in the expression for the single-site entanglement, 
and so the decay of the single-site entanglement is faster 
than linear in the region near A > 0. 

If there exist universal quantities related to the en- 
tanglement in critical quantum systems, then it is likely 
that they are derived from entanglement measures that 
satisfy additional properties beyond the set usually re- 
garded as 'essential' for an entanglement measure (see, 
for example, 1 40 , |43[). There are two main reasons why 



we make this assertion. The first arises from the inabil- 
ity of the single-site entanglement to distinguish between 
neighbouring and distributed entanglement. In order to 
distinguish between these differing scenarios, a good en- 
tanglement measure for critical quantum systems should 
take account of how the entanglement is shared out. The 
second reason is that, as we argue below, the single- 
site entanglement is not rescaleable. If a quantity is to 
be renormalizable it is necessary that it be rescaleable. 
That is, it must be possible to collect degrees of free- 
dom together, calculate the collective value of the quan- 
tity, and then rescale (or 'renormalize') the collective 
value. A renormalizable entanglement measure should 
be rescaleable in this way. 

We should be a little more precise in our definition of 
rescalcability for entanglement measures. Say we wish 
to calculate the bulk entanglement of a block of spins 
Si, «2) • • • j s m in a lattice with the rest of the lattice, L. 
If the entanglement measure §f (for example, could 
be the entanglement of formation) used to calculate this 
entanglement is to be rescaleable then, in the very least, 
it must satisfy the extensivity relation 



This inequality expresses the idea that the entanglement 
of a collection of spins with the rest of the lattice should 
be at least as great as the sum of the entanglements of 
each spin with L. If an entanglement measure does not 
satisfy the extensivity relation Eq. (5q) then it is not 
clear how to rescale the bulk value of the entanglement. 

Summarizing, the failure of the single-site entangle- 
ment to be universal may be due to the facts that: (a) 
it does not distinguish localized from distributed entan- 
glement; and (b) it is not rescaleable, in a sense that 
we can now make explicit. To do this, note first that it 
has previously been shown that the entanglement of for- 
mation does not satisfy Eq. ([38) Q. If we regard the 
single-site entanglement S as the entanglement of for- 
mation S = J^"(si,L) between a single spin s\ and the 
rest of the lattice L it seems unlikely that it will be a 
universal quantity. (There do exist other entanglement 
measures which reduce to the von Neumann entropy for 
pure states jjj [Tc| , fl3| , ffJJ . It is an open question whether 
they satisfy Eq. (p8|).) 

There are indications J48| , however, that the square of 
the concurrence is extensive. Perhaps a suitable gener- 
alisation of the concurrence will turn out to be the best 
quantity for studying universal properties of entangle- 
ment. Evidence that this is the case has recently been 
obtained by Osterloh et. al. Ej| where they found that 
a quantity related to the concurrence is universal for the 
transverse Ising and XY models. It would be interesting 
to investigate this behaviour and see if it arises because 
of the possible extensivity properties of the concurrence. 
Note, incidentally, that universal behaviour in the con- 
currence does not necessarily imply universal behaviour 
for the entanglement of formation, for the latter is only a 
function of the former in the special case of a two-qubit 
system. 

The determination of what entanglement is shared by 
two sites in the lattice requires a measure of the two-party 
entanglement present in mixed states. We will henceforth 
use the concurrence C to measure the two-party mixed- 
state entanglement between two spins. The concurrence 
of two spin-i particles may be calculated from their den- 
sity matrix g via the formula 



C{g) = max[0, Xi — A 2 — A3 — A4] 



(39) 



(si,s 2l ... ,s m :L)> &{sx : L) 



L). 
(38) 



where the A^ are the eigenvalues in decreasing order, of 
the Hermitian matrix R = */ 1 ^fgg^fg, and g = [a v ® 
a y )g*(a v ®a v ). The concurrence varies from C = for a 
separable state to C = 1 for a maximally entangled state 
& 

The two-site density matrices for the ground state of 
the XY model arc difficult to calculate when there is 
ground state degeneracy. This is because the magneti- 
sation (<7 X ) becomes nonzero as the phase- flip symme- 
try is broken, and it becomes necessary to include the 
correlation function (ctqCt^) in the operator expansion 
Eq. (18). The (0^0-*) correlation function is nonlocal 



when expressed in terms of the Jordan- Wigncr fcrmionic 
operators and there is no simple way to derive it from 
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FIG. 2: Nearest-neighbour concurrence C at zero temperature 

for the transverse Ising model FIG. 3: Next-nearest-neighbour concurrence C at zero tem- 

perature for the transverse Ising model 



other correlators. As a result of this difficulty we do not 
calculate the two-site density matrix for the ground state, 
instead, all two-site calculations are performed with re- 
spect to the thermal ground state. However, because the 
thermal ground state for the transverse Ising model takes 
the special form Eq. (p2[), it is possible to place bounds 
on the entanglement that can occur between two sites in 
a degenerate ground state. 

The entanglement between pairs of sites for the ther- 
mal ground state of the transverse Ising model shares 
many of the same features of the single-site entanglement. 
The entanglement, as measured by the concurrence, be- 
tween neighbouring sites and next-nearest neighbouring 
sites is shown in Fig. |2| and Fig. || respectively All other 
pairs have zero two-party entanglement because the cor- 
relation functions drop below the threshold for a positive 
concurrence. In both cases the entanglement rises from 
zero in the limits A = and A — > oo to a maximum value 
near the critical point A = 1. When A < 1 the ground 
state coincides with the thermal ground state so that the 
two-site entanglement results are the same in this case. 
Note that the maximum docs not occur exactly at the 
critical point A = 1. At first site this may appear to 
contradict our earlier conjecture that we expect entan- 
glement to be the g reates t at the critical point. In fact, 
as explained in Sec. Ill C , the reason for this is that the 
results here are for two-site entanglement, and are not 
inconsistent with the conjecture that the total entangle- 
ment in the lattice is a maximum at the critical point. 

It is interesting to see what effect the ground state de- 
generacy has on the two-site entanglement in the ground 
state. As mentioned, it is not possible to study the two- 
site entanglement for A > 1. Despite this difficulty, for A 
above the critical value, we can place a lower bound on 
the two-site entanglement in a degenerate ground state. 



This may be achieved by observing that the concurrence 
measure C is convex ltl5|L which means that 



(40) 



i=i 



i=l 



where pi is any probability distribution and pi a set 
of two-site density matrices. If we apply this inequal- 
ity to the thermal ground state, Eq. (|32|), we obtain 
C(por) < |C(tr o ^(|0+)(0+|)) + |C(tr (fr (|0-)(0-|)). The 
global phase flip is a local unitary operation, so that the 
concurrence of each term in the RHS of the inequality is 
the same, that is 



C(por) < C(tr fr (|0+)(0+|)). 



(41) 



In this way we see that the two-party entanglement in 
the ground state is at least as large as the two-party 
entanglement in the thermal ground state. 



C. Critical quantum systems and the constraints of 
shared entanglement 

The maximum value of the concurrence between neigh- 
bouring sites does not occur at the critical point. This 
seemingly contradicts the idea that the strength of the 
correlations is proportional to the entanglement, and that 
therefore the entanglement should be maximal at the crit- 
ical point. However, as we will discuss in this subsection, 
there are reasons based on the properties of shared en- 
tanglement to expect that this maximum should occur 
away from the critical point. 

It is well known that there are limitations to 
the amount of entanglement that may be distributed 
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amongst three or more subsystems [£5], Eq, [30| pi], |5^, 
Q. This class of problem, that is, the determination 
of how much two-party entanglement can be distributed 
amongst a given number of parties, is known as an en- 
tanglement sharing problem. The simplest example of 
this is the situation of three parties A, B and C. If A is 
maximally entangled with B then it is not possible for A 
and C or B and C to share any two-party entanglement. 
Entanglement sharing is relevant to the quantum phase 
transition in the transverse Ising model as it provides a 
fundamental bound on the amount of entanglement that 
may be distributed amongst the sites. The existence of 
such a bound means that as the overall entanglement in 
the lattice is increased, some sites become pairwise more 
disentangled. An example where this occurs is in a sys- 
tem approaching a critical point. 

As the critical point is approached in the transverse 
Ising model the correlation length begins to increase. 
What occurs physically is that each site develops entan- 
glement with its neighbouring sites. When the system 
gets closer to the critical point each site begins to develop 
entanglement with its next-nearest neighbours and so on. 
When the system is not at the critical point the entan- 
glement between a single site and the rest of the lattice 
is localised within some region because the correlations 
are exponentially damped for large enough separation. 
At the critical point this is no longer the case; there are 
appreciable correlations between a single site and every 
other site. However, the entanglement associated with 
this correlation must be distributed in such a way that 
it satisfies the constraints of entanglement sharing. We 
conjecture that the ground state at the critical point ac- 
tually saturates the bounds of entanglement sharing, so 
that it is maximally entangled in this sense. If this conjec- 
ture is correct, this would explain why the entanglement 
between neighbouring sites is not maximum at the crit- 
ical point. Initially, as A is increased, the entanglement 
between neighbouring sites increases first. When the sys- 
tem reaches criticality the entanglement is distributed 
to more remote pairs. If the ground state saturates the 
bounds of entanglement sharing this would have to occur 
at the expense of the two-party entanglement previously 
established between pairs of sites that are close. 

In the light of this interpretation it is interesting to 
compare the entanglement calculations for the transverse 
Ising model at criticality to the lattice calculations of 
Wootters and O'Connor [^5[ |mJ. In the critical case A c = 
1 the correlation functions for the transverse Ising model 
are known explicitly as functions of r [p8[ 
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FIG. 4: Nearest-neighbour concurrence C at zero temperature 
for the XY model 



where H(r) = \ r - 1 2 r ~ 2 . . . (r — 1). The concurrence at 
the critical point is nonzero for both r = 1 and r = 2 
where it is given by, respectively, 0.1946 and 0.0044. 
These values should be compared with the values ob- 
tained by O'Connor and Wootters in their study |30| of 
the concurrence in chains and rings of qubits. They max- 
imised the entanglement between nearest neighbours of 
a translationally invariant ring of spin-^ degrees of free- 
dom. Wootters and O'Connor were attempting to sat- 
urate the bounds of entanglement sharing by maximis- 
ing the entanglement of nearest neighbours subject to 
the symmetry of translational invariance. They found a 
maximal nearest-neighbour concurrence value of 0.4345 
for an infinite ring, which is greater than the critical value 
for the transverse Ising model. This result alone does not 
imply that the critical transverse Ising model is less en- 
tangled than the ring considered in |2^, Q , indeed, if the 
conjecture made in the previous paragraph is true then 
the ring would be much less entangled than the critical 
transverse Ising model. The reasoning for this is that 
the critical transverse Ising model is conjectured to max- 
imise the entanglement between all pairs subject to trans- 
lational invariance while the chains and rings of Woot- 
ters and O'Connor only maximise entanglement between 
nearest neighbours. One means of determining whether 
this is the case would be to calculate the correlation func- 
tion for the ring. On the basis of the arguments made 
in this study, we expect that the correlations will decay 
exponentially with separation for the ring. 

The entanglement in the thermal ground state of the 
general XY model may be calculated simply, follow- 
ing the method outlined in Sec. II. Following Barouch 
fS7f , which is where the correlation functions Eq. (|2~j]), 
Eq. (|22]), and Eq. ( p3| ) were calculated, only the re- 
gion < 7 < 1 is considered here. The concurrence 
between nearest-neighbour and next-nearest neighbour 
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FIG. 5: Next-nearest-neighbour concurrence C at zero tem- 
perature for the XY model 



sites is shown in Fig. g| and Fig. |5j respectively. The con- 
currences are a complicated function of the parameters, 
reflecting the competition between the various different 
noncommuting terms in the Hamiltonian as the parame- 
ters are varied. 

The completely isotropic limit, 7 = 0, is the most in- 
teresting parameter region besides the transverse Ising 
model. Direct calculation along the lines already pre- 
sented shows that two-party entanglement exists between 
all pairs for all separations at this point. Wootters ]55| ] 
has made a study of the correlations in one- and two- 
dimensional lattices and he has found interesting connec- 
tions between the two-party correlations in the isotropic 
XY model and the bounds of entanglement sharing. Fur- 
ther investigations along these lines could provide ev- 
idence that critical quantum lattice systems are maxi- 
mally entangled in the sense of entanglement sharing. 



IV. THERMAL ENTANGLEMENT IN THE 
TRANSVERSE ISING MODEL 

In this section we discuss the entanglement present in 
the thermal state of the transverse Ising model. We find 
that the largest amount of entanglement is present in 
the parameter region close to the critical point. This 
region is found to correspond with the quantum critical 
region introduced by Sachdev (|l6||, pg. 58). We also find 
parameter values for which the entanglement increases 
as the temperature is increased. Finally, we discuss the 
persistence of quantum effects in the thermal state as the 
temperature is increased. 

It is desirable to determine when a condensed-matter 
system will behave quantum-mechanically. This is par- 
ticularly important because the validity of various ansatz 
methods depends on whether they take account of possi- 



FIG. 6: Nearest-neighbour concurrence C at nonzero temper- 
ature for the transverse Ising model 




FIG. 7: Next-nearest-neighbour concurrence C at nonzero 
temperature for the transverse Ising model 



ble quantum effects. When a system is in its ground state 
quantum effects will certainly be important, as evidenced 
by the quantum phase transition in the XY model. The 
zero-temperature calculations of the last section repre- 
sent a highly idealised situation, however, and it is un- 
clear whether they have any relevance to the system at 
nonzero temperature. It turns out that the properties 
of a quantum system for low temperatures are strongly 
influenced by nearby (in parameter space) quantum crit- 
ical points |l6|, [l?]]. It is tempting to attribute the effect 
of nearby critical points to persistent mixed-state entan- 
glement in the thermal state. In order to investigate 
this, we calculate the two-party entanglement present at 
a nonzero temperature T. 
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The two-site density matrices constructed in Sec. II are 
valid for all temperatures. Using these matrices it is pos- 
sible to study the purely two-party entanglement present 
at thermal equilibrium because the concurrence measure 
of entanglement can be applied to arbitrary mixed states. 
The regions where there is appreciable two-party entan- 
glement give at least a partial indication of where quan- 
tum effects may be important. We again emphasise the 
transverse Ising model for this section. The influence 
the critical point has on the entanglement structure at 
nonzero temperatures is particularly clear for this model. 

The entanglement between nearest-neighbour and 
next-nearest neighbour sites in the Ising model at 
nonzero temperature appears in Fig. |^ and Fig. [?] respec- 
tively. The entanglement is nonzero only in a certain 
region in the fc^T — A plane. It is in this region that 
quantum effects are likely to dominate the behaviour of 
the system. The entanglement is largest in the vicinity 
of the critical point A = 1, ksT = 0. This region cor- 
responds, approximately, to the quantum critical regime 
identified by Sachdcv Sachdev found, by using a 

very different argument, that quantum effects would be 
important in this regime. The correspondence of these 
two regions provides evidence that the entanglement con- 
tent plays an important role in the emergence of quantum 
behaviour in naturally occurring quantum systems. 

There are two notable features of the two-site thermal 
entanglement results. The first feature is that, for cer- 
tain values of A, the two-site entanglement can increase 
as the temperature is increased (eg. A = 1.4, Fig. ^|). This 
effect has previously been observed in finite-size calcula- 
tions JlJ, g(| for the Heisenberg model. The occurrence 
here of the same effect implies that it is not an artifact 
of the truncation of a lattice. The second feature is the 
existence of appreciable entanglement in the system for 
temperatures ksT above the ground state energy gap 
A. It has been argued |l7| that quantum systems behave 
classically when the temperature exceeds all relevant fre- 
quencies. For the transverse Ising model the only rele- 
vant frequency is given by the ground state energy gap 
A = huj. The presence of entanglement in the system for 
temperatures above the energy gap indicates that quan- 
tum effects may persist past the point where they are 
usually expected to disappear. 

A comparison should be made between the results ob- 
tained here and the numerical calculations of concur- 
rence in the Ising model on a finite number of sites p8| . 
The calculations that were performed in ps|| were im- 
plemented on a maximum of 7 sites. The concurrence 
between nearest neighbours obtained by Gunlycke et al. 
(Figs. 2 and 5 of p8| ) is in qualitative agreement with 
the results obtained here. However, as there is no phase 
transition for the finite size Ising model the dominance 
of the critical point was not as sharp in the calculations 
of©- 



V. SUMMARY AND FUTURE DIRECTIONS 

The one- and two-party entanglement present in the 
ground and thermal states of the XY model has been 
calculated. It should be stressed that the calculations in 
this study are analytic and, furthermore, they are for the 
thermodynamic limit of a quantum lattice system. 

We have argued that the critical point of a quantum 
lattice system corresponds to the situation where the lat- 
tice is maximally entangled. Evidence for this conjec- 
ture was found in the single-site entanglement results for 
the ground state of the transverse Ising model. We have 
also argued that the constraints of shared entanglement 
are important for critical quantum systems, and we have 
found possible evidence of such constraints playing a role 
in the two-party entanglement results for the transverse 
Ising model. The entanglement present at thermal equi- 
librium was also studied, and an approximate correspon- 
dence between the quantum critical regime identified by 
Sachdev and the regions where the two-party entangle- 
ment is nonzero was found. Parameter values where the 
entanglement increases as the temperature is increased 
were also found. 

We have focused on the transverse Ising model 
throughout this study, although the calculations pre- 
sented also cover the XY model. The transverse Ising 
model is interesting because it is the simplest system to 
exhibit a quantum phase transition, and it is relatively 
easy to identify the structure of the entanglement present 
in this system. The importance of the critical point in 
this system is also particularly clear. The XY model has 
many parameter regimes where it behaves differently, so 
it is very likely that more interesting phenomena may be 
found in other parameter regions. 

Entanglement calculations in this study have been re- 
stricted to time-independent scenarios. However, the dy- 
namic correlation functions have been calculated for the 
Ising and XY models for certain values of A. It is possible 
and may be interesting to calculate the time evolution of 
the entanglement in these models and thus identify truly 
quantum dynamics. 

The calculations in this study arc intended as a point of 
reference for the development of an understanding of the 
entanglement in critical quantum systems. Rather frus- 
tratingly, the present incomplete understanding of en- 
tanglement measures has prevented us from performing 
many of the calculations we would like to do in order to 
check the many conjectures made in this paper. Further 
progress on the general quantitative theory of entangle- 
ment should enable these conjectures to be checked in the 
future. We believe that entanglement plays a central role 
in the emergence of long-range correlations at the critical 
point of such systems, and that a fruitful interplay be- 
tween the theory of entanglement and critical quantum 
phenomena may result from further study. In particu- 
lar, it would be interesting to make universal statements 
about the character of entanglement at the critical point, 
and to examine whether the constraints of entanglement 



13 



sharing impose physical limitations on the behaviour that 
can occur in such a system. 

Note added: As this paper was nearing completion we 
learnt of related work done independently by Osterloh et. 
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